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The mean inetallicity of the thick disk of the Galaxy is 0.5 dex higher than that of 



^ ! Abstract 

o 

^\| ■ the damped Lya systems. This has been interpreted to argue that stars in the former 

do not arise out of gas in the latter. Using new metallicity and H I column-density data 

we show the metal-rich damped systems do contain sufficient baryons at the thick-disk 

metallicity to account for the stellar masses of thick disks. Comparing our kinematic 

data with the metallicities we show that damped Lya systems exhibiting the largest 

_ profile velocity widths, Av, span a narrow range of high metallicities, while systems with 

^ ■ small At; span a wider range of metallicities. This is naturally explained by passage of 

(^ . the damped Lya sightlines through rapidly rotating disks with negative radial gradients 

in metallicity. The systematically lower A'^(H I) of systems with high Av indicates (a) the 

gaseous disks have centrally located holes, and (b) an apparent inconsistency with the 

^ . protogalactic clump model for damped Lya systems. The higher metallicity of systems 

H I with low A^(H I) further implies that stars rather than gas dominate the baryonic content 

of the most metal-rich damped systems. 
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1 INTRODUCTION 

The collapse of a spheroidal protogalaxy to the centrifugally supported disk of the Galaxy was 
inferred from correlations between the metallicities and kinematics of old stars in the solar 
neighborhood (Eggen, Lynden-Bell, and Sandage 1962). But subsequent studies have not 
sorted out the sequence of events leading to the formation of stellar populations comprising 
the halo, the thick disk, and the thin disk (e.g., Majewski 1993). The damped Lya absorption 
systems, a population of H I layers widely believed to be the gaseous progenitors of current 
galaxies (see Wolfe 1995), provide an independent perspective for studying these events 
because (a) they occur in objects comprising the bulk of the galaxy population at high 
redshifts, and (b) the ranking of redshifts yields an unambiguous time sequence. Detected in 
the redshift interval z = [0,4.5] the damped Lya systems trace the evolution of neutral gas 
in galaxies from their protogalactic phase to the present. However, the chemical properties 
of the gas may be incompatible with those of existing stellar populations. At z = [1.6,4.5] 
the metallicity of the gas is low compared to the thin disk metallicity (Pettini et al. 1997) 
indicating that stars in the thin disk do not arise directly from high-z damped systems 
(Lanzetta et al. 1995). The possible enhancement of alpha-rich elements suggests the gas 
gives rise to halo stars (Lu et al. 1996), but the kinematics of the gas are inconsistent with 
this hypothesis (Prochaska & Wolfe 1997). More recently Pettini et al. (1997) argued that 
the metallicities of the damped systems are too low to explain the thick disk (Gilmore et al. 
1989; Carney et al. 1996). 

In this letter we reconsider the scenario in which star formation in damped Lya systems 
results in the formation of the thick disk. Combining metallicities and column densities with 
new kinematic data obtained with the Keck I 10 m telescope we suggest a plausible scenario 
in which the thick disk forms out of damped Lya gas. 

2 COSMIC METALLICITY DEPENDENCE ON 
BARYON DENSITY 

We wish to find whether the mass content and metal abundances of gas in damped Lya 
systems can account for the mass density and metallicities of thick stellar disks. Define the 
cosmic metallicity < Z > = ^metais/^g where ^metais and Qg are the comoving densities of 
metals and neutral gas in damped Lya systems (Lanzetta et al. 1995). Let the number of 
damped systems in the metallicity and column-density intervals {Z',Z'+dZ') and (A^,A^+dA^) 
be given by h{Z', N)dZ'dN. The latter is related to the frequency distribution of H I column 
densities by f{N)=J h{Z',N)dZ' (Lanzetta et al. 1995), and the frequency distribution of 
metallicities by g{Z')=J h{Z',N)dN. Suppose h{Z',N) spans the metallicity interval Z' = 
[Zmin, Zmax] and columu-density interval N = [Nmin, Nmax], and ^g{Z) is the density of 
damped Lya baryons in the metal-rich subinterval Z' = [Z,Zmax]- Then ^g{Z) and the 
corresponding < Z[Z) > are given by 
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Zi)6{N — Ni), where the sum extends over all the iVj, Z'^ pairs in the sample. As a result 

Qk = Qg X ^ , <Z^>= ^=i- , (3) 
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where the indices i=l, k, and imin correspond to Zmax, Z, and Zmin- Because the sums in eq. 
(3) are over an array of damped Lya gas layers ordered according to decreasing metallicity, 
< Z{Z) > decreases with decreasing Z while ^g{Z) increases. We can determine ^g{Z) 
corresponding to the mean metallicity of the thick disk, provided the latter is less than Z^ax- 

To determine < 2{Z) > as a function of ^g{Z) we turn to the [Zn/H], A^(H I) pairs 
that Pettini et al. (1997) acquired for 34 damped systems, where Z=Zq10'^'^'^'. We focus 
on Zn rather than Fe as a metallicity indicator, because this is the largest recorded sample 
of damped Lya metal abundances, and Zn is less depleted than Fe by dust which may be 
present (Fall & Pei 1995). We select the 27 pairs in the redshift range z =[1.6,3.0]. Systems 
with z > 3.0 are excluded, since the metallicities in this redshift range are systematically lower 
than those of the thick disk. The sample comprises 16 systems with detected [Zn/H] and 11 
with upper limits. Two of the detections, which come from our Keck HIRES observations, 
replace the upper limits of Pettini et al. (1997). 

We used equation (3) to determine the points in Figure 1 which plots [<Zn/H>] (= 
\og{Zk/ZQ)) vs Vtk for k = [l,imin]- We let Qg = 0.003, the value inferred by Storrie-Lombardi 
& Wolfe (1997) in the redshift interval z = [1.8,3.5] for Hq = 50 km s~^ Mpc~^ (which is 
adopted throughout this paper). The circles were computed by letting the upper limits equal 
the true values of [Zn/H]. In this case ilg{Z) = 0.0004 when [<Zn/H>] = — 0.6, the mean 
metallicity of the thick disk (Carney et al. 1996). The triangles and squares were computed 
by equating the upper limits minus 0.5 and 1.0 with the true values of [Zn/H]. In both cases 
^g{Z) ^ 0.0003 when [<Zn/H>] = — 0.6, indicating the result is robust. Assuming bulges 
and disks contribute equally to the density of visible matter (Schechter & Dressier 1987), 




Figure 1: Log of mean cosmic metallicity vs the comoving density of damped Lja baryons, 
for damped systems with metalhcities Z = [Z, Z^ax]- Circles computed assuming upper hmits 
equal true values of [Zn/H]. Triangles and squares computed assuming upper limits minus 0.5 
and 1.0 equal true values of [Zn/H]. In latter cases all baryons are assumed to be gas. Thus 
eq. (3) is used to compute Qdlab{Z) = Q/. and [<Zn/H>] = log(< Z^ >/Zq). Stars include 
correction for the presence of baryons in stars. In this case eqs. (4) —(6) are used to compute 
^dlab{Z) and [<Zn/H>]. Vertical and horizontal dashed lines correspond to cosmic density 
and metallicity of the thick stellar disk. 

which is given by 0.0054 (Gnedin & Ostriker 1992), and that the mass of the thick disk is 
0.1 times that of the thin disk (Majewski 1993) we find the thick disk mass density, ^thick = 
0.00027. Although the error bars associated with QtMck are of order 50%, it is reasonable to 
conclude that the damped Lja systems contain sufficient baryons to account for the masses 
of thick stellar disks (see Figure 1). 



3 KINEMATICS, METALLICITIES, AND STARS 

To learn more about the metal-rich damped systems we turn to the kinematics of the gas. 
Analysis of the velocity profiles of weak metal lines in over 30 damped Lya systems shows 
the frequency distribution of profile velocity widths, Av, and other statistics that test for 
asymmetries exhibited by the profiles are consistent with absorption by thick disks with 
rotation speeds Vrot ~ 250 km s^^(Prochaska & Wolfe 1997). The CDM simulation of Haehnelt 
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Figure 2: (a) metallicity vs Av for damped systems with z = [1.6,3.0]. Filled and empty 
circles correspond to detections and upper limits, (b) logA^(H I) vs Av. 



et al. (1997), in which infall, random motions, and rotation of protogalactic clumps contribute 
equally to Av, may be likewise consistent. Here we focus on rotating disks. 

Figure 2a plots 17 [Zn/H], Av pairs drawn from our kinematic sample with z =[1.6,3.0]. 
The figure shows that systems with high Av and low metallicity are not detected in this 
redshift range. Specifically metallicities [Zn/H] < —1.0 are absent in all 5 systems with Av 
> 120 km s~^ , but present in 7 out of 12 systems with Av < 120 km s~^ . The effect is 
real and is not an artifact due to observational selection, since systems with high Av and low 
metallicity are detected at z > 3.0. Nor is dust likely to be a contributing factor, since dust 
would remove metal-rich rather than the metal-poor systems missing from Figure 2a. The 
reality of this effect is further supported by its presence in a [Fe/H] vs Av diagram. We also 
find possible evidence for a correlation between [Zn/H] and Av exceeding 3.5o" significance 
when the true [Zn/H] equal the upper limits minus 1.0. 

The systematic pattern in Figure 2a can be explained by negative radial gradients in 
metallicity. Monte Carlo simulations of absorption profiles produced by sightlines penetrating 
randomly oriented disks indicate that a necessary condition for large Av is a small impact 




2 

b 



Figure 3: Distribution of impact parameters resulting from the numerical simulation described 
in text. Solid curve corresponds to impacts leading to Af > 120 km s^^ , and dotted curve 
to Af < 120 km s^^ . Impact parameters in units of radial scale-length, R^. 



parameter. This is evident in Figure 3 which plots the distribution of impact parameters, b 
(where b is in units of radial scale length, Rd, of an assumed exponential gas distribution), 
resulting from simulating identical exponential disks with rotation speed Vrot = 250 km s~^ 
and vertical scale-height h = O.SRd- Whereas 86% of impacts leading to Av > 120 km s~^ 
are confined to 6 < 1, 1% are at 6 > 2. Therefore, the absence of low metallicities at Av 
> 120 km s~^ requires high element abundances at small radii. On the other hand 26% of 
impacts leading to Av < 120 km s~^ are at 6 < 1, while 16% are at 6 > 2. The wide range of 
impact parameters can explain the broader distribution of metallicities at Av < 120 km s~^ , 
if impacts at large b yield low metallicities, i.e., if metallicity decreases with radius. None of 
these results changes significantly when we use a more realistic model in which Vrot is drawn 
from a distribution of rotation speeds characterizing present-day spiral galaxies. 

Damped Lya systems with large Av also exhibit systematically lower A^(H I). This is 
shown in Figure 2b which plots logA^(H I) vs Av for 29 damped Lya systems drawn from our 
kinematic sample. Whereas 1 out of 10 systems with Av > 120 km s'^ has logA^(H I) > 20.6, 
12 out of 19 systems with Av < 120 km s~^ have logA^(H I) > 20.6. Figure 2b includes systems 
with z > 3.0, because the effect is independent of redshift. Suppose the gas distribution has 
a central hole. At high Av the impact parameters are so small that the sightlines encounter 
the low column densities present at small radii. A wider range of iV(H I) occur at low Av, 



because the sightlines sample a broader range of impact parameters. Preliminary results from 
simulations with central holes are in better agreement with the logA^(H I) vs Af data than 
standard exponential disks. 

Deficiency of neutral gas occurs often in the central regions of spiral galaxies (Broeils & 
van Woerden 1994), the same regions where enhancements in metallicity are also common 
(Edmunds & Pagel 1984). Thus, there is empirical support for the idea that damped Lja 
systems comprise gaseous disks with central holes and negative radial gradients in metallicity, 
if they evolve into current spirals. The increased metallicity is a signature of enhanced star 
formation which also helps to explain the deficit of gas, either through direct gas consumption 
or the loss of gas through energetic outflows from supernovae. In either case a significant 
fraction of baryons in the metal-rich damped systems may be locked up in stars. As a result 
the expression for cosmic metallicity in eq. (3) will underestimate the contribution from the 
gas-poor metal-rich systems. 

When stars are present, Qk and < Zk > are given by 

n, = n^x^^i- , <Z, >=i=^- , (4) 

3=1 i=i 

where Vt^ is the comoving mass density of baryons in stars plus gas. Because N^ and Z^ are 
the column density and metallicity of matter in stars, < Zk> is the comoving mass density of 
metals in stars plus gas divided by Vtk- Although < Zk> vn eq. (4) differs from the standard 
definition for metallicity, it is the appropriate quantity, because metals in stars as well as gas 
in damped Lya systems supply metals to stars comprising the current thick disk. And the 
thick disk metallicity is inferred solely from stars. 

To solve eq. (4) we first adopt the chemical evolution model of Larson (1972) to compute 
the fraction of baryons in stars. The model assumes the star formation rate is balanced by 
the rate of mass infall to the disk, and as a result the gas content of the galaxy does not 
change. This agrees with the observed constancy of ^g{z) in the redshift range z = [1.6,3.3] 
(Storrie-Lombardi & Wolfe 1997) in which stars of the thick disk are assumed to form. We 
have 



N:/N, = In 



y + Zf-Z, 



init 



y + Zf-zi 



(5) 



where y is the chemical yield, and Zf and Zinu are the metallicities of the infalling material 
and of the 'initial' disk aX z > 3.0. We determine Zf from the constraint 



ZtxNt + {Z[- Z,^,t)^N,=yxNt , (6) 

(see Tinsley 1980). We combined eqs. (4) — (6) to compute \og{Zk/ZQ) versus Vtk in the 
presence of stars. The solution, shown as stars in Figure 1, was computed assuming y = O.BZq, 



Zinit = Zf = O.OIZq. We estimated Zinn and Zf from the lowest metallicities found for 
damped Lya systems ai z > 3.0, and y from standard models (Tinsley 1980). In this case 
Vlk = 0.0008, i.e., SQthick, when log{Zjs/ZQ) = — 0.6. The increase in Qj^ at the metallicity 
of the thick disk results from the significant stellar corrections in baryonic mass for the metal 
rich damped Lya systems. 

4 DISCUSSION AND CONCLUSIONS 

Our results suggest that contrary to previous claims (Pettini et al. 1997) the damped Lya 
systems contain more than enough baryons at suitable metallicity and rotation speed to form 
thick stellar disks in spiral galaxies. The kinematic/metallicity data further imply (a) stars in 
the thick disk form in the inner metal-rich regions of rapidly rotating gaseous disks, and (b) 
these stars may dominate the baryonic content of the most metal rich damped systems. We 
conjecture that vertical contraction of the metal-rich, thick gaseous disk leads to the formation 
of the inner thin disk. The metal-poor gas of the outer disk, i.e., the ~ 90 % of the damped 
Lya baryons that remain in gas after the formation of the thick disk, could supply the 
remaining thin-disk mass through radial contraction driven by angular momentum transport 
mediated by high-amplitude spiral density waves (Roberts & Shu 1972). The metallicity of 
this gas may increase as a result of mass loss by stars in the inner regions. 

By contrast stars forming in the protogalactic clumps considered by Haehnelt et al. (1997) 
end up in bulges and halos rather than rotationally supported structures, because the motions 
of the merging clumps are not dominated by rotation. Formation of the thick disk occurs at 
z < 1, after merging ceases, and when the mean metallicity equals —0.6. But the age of the 
thick disk is unlikely to be less than 12 Gyr (Carney 1997) which exceeds the lookback time 
to 2; = 1 in all A = cosmologies and spatially flat cosmologies in which ^a < 0.8, when Hq 
> 50 km s~^ Mpc~^. Edge-leading asymmetries of the velocity profiles arise in this model 
because clumps with the highest volume and column density move fastest with respect to the 
surrounding gas. This may occur because ram pressure deceleration by ambient gas is less 
effective in decelerating denser clumps. The predicted correlation between A^(H I) and Av is 
in conflict with the detection of systems having high Av (i.e., > 120 km s~^ ) and low logA^(H 
I) (< 20.6) (see Figure 2b). 
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